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Orbifold equivalence and the sign problem at finite baryon density 
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We point out that SO(2N c ) gauge theory with Nf fundamental Dirac fermions does not have a 
sign problem at finite baryon number chemical potential fis ■ One can thus use lattice Monte Carlo 
simulations to study this theory at finite density. The absence of a sign problem in the SO(2N c ) 
theory is particularly interesting because a wide class of observables in the SO(2N c ) theory coincide 
with observables in QCD in the large N c limit, as we show using the technique of large N c orbifold 
equivalence. We argue that the orbifold equivalence between the two theories continues to hold at 
finite \ib provided one adds appropriate deformation terms to the SO(2N c ) theory. This opens up 
the prospect of learning about QCD at finite fiB using lattice studies of the SO(2N c ) theory. 



The properties of QCD at high baryon densities have 
long been a subject of intense interest. Apart from its 
intrinsic theoretical appeal, this subject is important in 
astrophysics, especially in the study of neutron stars. Be- 
cause of asymptotic freedom, the behavior of QCD at 
asymptotically high chemical potential for baryon num- 
ber hb is well understood theoretically, and QCD be- 
comes a color superconductor as fis — > oo [lj. At more 
phenomenologically realistic densities, QCD is strongly 
coupled, and thus not amenable to controlled analytic 
treatment. Lattice Monte Carlo simulation is very use- 
ful at fj>B = 0. However, it runs into trouble at [Ib 7^ 
due to the fermion sign problem: the fermion determi- 
nant becomes complex, rendering importance sampling 
exponentially difficult. 

Over the years, several gauge theories that do not suffer 
from sign problems at finite density have been explored. 
The main examples are QCD with an isospin chemical 
potentially, two-color QCD0 g], and adjoint QCDg]- 
6 . However, while interesting, these theories have many 
qualitative differences from N e = 3 QCD, such as e.g. 
explicitly broken flavor symmetry in the first case. 

Here, we propose a path to study QCD at [Ib 7^ in 
the large N c limit using lattice Monte Carlo methods. 
Large iV c QCD [7] gives many insights into nonpertur- 
bative strong interactions at zero \xb'- it is often a good 
approximation to our N c — 3 world. The extent to which 
the finite-density large N c world is a good approximation 
to the N c = 3 world is discussed in, for instance, Refs. [8]. 

Our proposal rests on two observations. The first is 
that all representations of SO(2N c ) are real, and as a re- 
sult SO(2N c ) gauge theory with Nf fundamental Dirac 
fermions does not have a sign problem at /i^ 7^ 0. This 
alone already makes the theory worth studying, espe- 
cially because the SO(2N c ) theory shares a number of 
qualitative features with N c > 3 QCD. For example, it 
has 2A^ c -valence-quark baryons [9]. 

Our second observation is that the connection between 
the SO(2N c ) theory and SU{N C ) QCD is in fact quan- 
titative. We show that SU(N C ) gauge theory with Nf 



fundamental Dirac fermions (i.e., large N c QCD) can be 
obtained as an orbifold projection of the SO(2N c ) theory. 
Large N c orbifold equivalence [T0HT2] then guarantees 
that all correlation functions of operators in the "neu- 
tral" sector (i.e., invariant under the symmetry used for 
the projection) coincide in both theories to leading order 
in the 1/N C expansion, provided the symmetries used in 
the projection are not spontaneously broken. The nec- 
essary symmetries are unbroken at fiB = 0, and thus 
the SO(2N c ) theory and large N c QCD have coincid- 
ing correlation functions for a broad class of operators. 
The equivalence should continue to hold at //g 7^ 0, pro- 
vided one adds certain deformation terms to the SO(2N c ) 
theory which protect the orbifold symmetry, but do not 
otherwise affect the connection of the theory to large N c 
QCD. We show that there exist deformations that pro- 
tect the orbifold symmetry at least for *C ^qcd, and 
likely for larger /j,b as well, all while keeping the the- 
ory sign problem free in the chiral limit. The existence 
of a sign-problem-free theory equivalent to finite-density 
large N c QCD is unexpected and quite remarkable. 

The orbifold equivalence between the two theories at 
/j,b = can be checked using lattice simulations, as can 
the question of whether the necessary symmetries are 
protected at large \ib- If the proposal passes these checks, 
the SO(2N c ) gauge theory may be used to perform non- 
perturbative studies of large N c QCD at finite density. 

SO(2N c ) gauge theory. The SO(2N c ) 4D gauge the- 
ory with Nf flavors of Dirac fermions (in Euclidean sig- 
nature) is 

Nf 

Cso = -A- Tr F^+Y^qai^D^+m. + ^Bl^qa (1) 

QSO a=l 

where F^ v is the SO(2N c ) field strength, Z? M = d^ + iA^, 
q a is a Dirac fermion in the fundamental representation 
of SO(2N c ), and m q and hb are the quark mass and 
chemical potential. = A^U, where the ti are the 

generators of SO(2N c ); we take Tr trfj — Sij. 

When m q = Hb = 0, Eq. Q has an SU(N f ) L x 



2 



and Nf 



SU(Nf) R x U(1)b X U(\)a chiral symmetry at the clas- 
sical level, just like SU(N C ) QCD. However, the chi- 
ral symmetry of the theory is actually larger than this, 
since SO(2N c ) is a real gauge group; classically it ex- 
tends to U(2N f ) [13] . U(1) A C U(2N f ) is anoma- 
lous as usual at finite N c . The chiral condensate qq 
breaks to SO{2N f ) D SU(N f ) v . The resulting Nambu- 
Goldstone bosons (NGBs), with mass ~ y'm,, live 
on SU(2N f )/SO(2N f ). Some of the NGBs, the "pions", 
are pseudoscalars that couple to g a 75<7fa, while the oth- 
ers, which we will refer to as baryonic pions, are charged 
under U(1)b- The baryonic pions are parity even [14] 
and couple to color-singlet operators of the form S a b and 
S^ ab , where S a b — q^C^qb and C — 7472 is the charge 
conjugation matrix satisfying Cj^C^ 1 = — 7*. 

There are Nf(Nf + 1) baryonic pions in the theory, 
1 pions with no baryon number. The SO(2N c ) 
theory also contains baryon-number-charged cousins of 
other mesons normally encountered in QCD. 

Now consider turning on fis 7^ 0. Since in the chiral 
limit the baryonic pions are the lightest particles charged 
under U(\)b, once > m„/2, one would expect the 
system to undergo a second-order phase transition to a 
phase with a nonzero density of baryonic pions. In fact, 
on general grounds, one expects that the baryonic pi- 
ons will Bose condense. In other theories with "baryonic 
pions", namely 2-color QCD and adjoint QCD, explicit 
chiral perturbation theory (xPT) calculations show that 
this does indeed happen @]. We expect the same in 
the SO(2N c ) theory [14]. The baryonic pion conden- 
sate breaks U{1)b — > ^2- The breaking of U(1)b at 
Ms > m Tr/2 in the SO(2N c ) gauge theory is in sharp 
contrast to the way SU(N C ) QCD behaves, where there 
are no baryonic pions to be condensed. We return to this 
crucial point below, in the context of orbifold projections. 

Orbifold projection to SU(N C ). To perform an orb- 
ifold projection, one identifies a discrete subgroup of the 
symmetry group of the "parent" theory, which for us 
is the SO(2N c ) theory, and sets to zero all of the de- 
grees of freedom in the parent theory that are not invari- 
ant under the discrete symmetry. This gives a "daugh- 
ter" theory, which in this case turns out to be large N c 
QCD. The orbifold projection uses a I2 subgroup of the 
SO(2N c ) x U(1) B symmetry of the SO{2N c ) theory. 

To define the orbifold projection, take J G SO(2N c ) 
to be given by J = io-% ® 1jv c ; In is an N x N identity 
matrix. [For earlier work on projections from SO(2N c ) 
to SU(N C ), see [T5].] J generates a Z 4 subgroup of 
SO(2iV c ). Next, let u = e"/ 2 € U(1) B generate a Z 4 
subgroup of U(1)b- The action of J and ui on A^, q a is 



Afj, y J A^J 



Qa 



-0JJq a , 



(2) 



A^ can be written in N c 
( A A 



N r blocks as 



(3) 



where fields with an 'A' ('S") superscript are antisym- 
metric (symmetric) matrices. Under the Z2 symmetry, 
A£,D^ are even while B^,C A are odd, so the orbifold 
projection sets B A = C„ = 0. So 



J\proj 



D s , A? 



If one defines a unitary matrix 



1 



f 1 ' 



~ V2 \ 1n c ~ il N c 



then 



PA projp- 



-Al 
A u 



(4) 



(5) 



(6) 



where A^ = + iA A is a U(N C ) gauge field. However, 
the difference between U(N C ) and SU(N C ) is a 1/N% cor- 
rection. The gauge part of the action of the orbifold- 
projected parent theory is thus simply 



C 



gauge, proj _ 



Tr T^T^. 



(7) 



where J-^ v is the SU(N C ) field strength. 

Now consider the effect of the orbifold on q a . Writing 
(A+, A~) T = (Pq a ) T , the action of the Z 2 symmetry is 
just (A+,A~) T — > (— A+,A~) T . The projection consists 
of setting A+ = 0. 

The action of the daughter theory is the action of the 
parent theory after the projection, with a rescaled cou- 
pling constant g S u = 9SO [E] 



1 



C=-^Tr J%+Y, ^ (-fV* + m i + ipa (8) 

^9su 



generating a Z2 subgroup of SO(2N c ) x U(1)b- 



where T^v is the field strength of the SU(N C ) gauge field 
A,t = Dl + iA£, i> a = A", and X> M = d M + iA^. This 
is an SU(N C ) gauge theory with Nf flavors of funda- 
mental Dirac fermions. So the orbifold projection relates 
SO(2N c ) gauge theory to large N c QCD. 

Neutral sector. The claim of orbifold equivalence 
is that the connected correlators of neutral operators 
in orbifold- equivalent parent and daughter theories will 
agree at large N c . We define neutral operators to be 
those that are invariant under the projection symmetry. 
Color-singlet gluonic operators in the SO(2N c ) theory 
are neutral, and are mapped to C-even gluonic operators 
in SU(N C ) theory by the projection. 

For fermionic observables, things are more subtle. 
Some examples of fermion bilinears that survive the pro- 
jection are given in Table I. An example of a bilinear that 
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SO(2N c ) theory 


SU{N C ) theory 


q a qb 


Tpalpb 




ipa-y^tpb 




ipa-y^d^b 



TABLE I. Examples of fermion bilinears whose correlation 
functions match between the two theories at large N c . 



does not survive the projection is q^C^qb, the baryonic 
pion operator; the same is true for all of the baryonic 
mesons. This is because all such operators have charge 
— 1 under the Z2 projection symmetry. Thus correla- 
tion functions involving baryonic mesons have no coun- 
terparts in large N c QCD. Note that this implies that 
the counting of neutral NG bosons matches in the two 
theories which is an important sanity check on the 
large N c equivalence. 

Validity of the equivalence. Given an orbifold projec- 
tion between two theories, the correlation functions of 
neutral operators will agree in perturbation theory [TU], 
For an orbifold equivalence to hold nonperturbatively, 
the symmetries used in the projection must not be spon- 
taneously broken [IT]. While the proofs of Refs. [UJ need 
be generalized to apply to projections involving funda- 
mental fermions|14l [Jj)][T7], it is natural to conjecture 
that the symmetry realization condition remains the key 
to the nonperturbative validity of the equivalence. 

In our case, U{1)b — > 7t>i when (is > "^7r/2 due to 
baryonic pion condensation. This breaks the projec- 
tion symmetry, destroying the equivalence. All is not 
lost, however. One can add a deformation term to the 
SO(2N c ) theory to prevent baryonic pion condensation 
and protect U(1)b- The simplest choice is 



so 



cf 

A2 



(9) 



where A ~ ^qcd- The deformed theory should be 
viewed as an effective field theory, defined with an im- 
plicit (lattice) cutoff, since the deformation is an irrele- 
vant operator. Large N c factorization implies that the 
deformed system would pay an energy cost 0(c 2 ) for the 
formation of a baryonic pion condensate. Thus provided 
c is large enough, the deformation prevents baryonic pion 
condensation, saving the validity of the orbifold equiva- 
lence between the deformed SO(2N c ) gauge theory and 
large N c QCD. 

By construction, deformation terms do not survive the 
orbifold projection to QCD, so the value of c does not 
affect the correlation functions of neutral operators so 
long as the equivalence holds. Much like the double- 
trace deformations used to prevent [18] center-symmetry 
breaking in Eguchi-Kawai reduction [15] , our deforma- 
tion terms hide themselves once they do their job. 

The sign problem. Consider the undeformed SO(2N c ) 



theory. Then the Dirac operator D = Jp + rn q + hbJa 
satisfies Cj 5 D(C'f)- 1 — D*. If the lattice form of 
the Dirac operator also has this symmetry, then if ip 
satisfies Dip = Xip, D^C^ip*) = A*(7 5 C"V), and 
tp, 7 5 C~ 1 (p* are orthogonal [S]. So eigenvalues form pairs 
(A, A*), and hence det(D) > 0, even at [Ib ^ 0. 

For simulations of the deformed theory, the action 
must be made quadratic in q, which can be arranged 
by 'integrating in' auxiliary fields. For general values of 
m ? ,/ig, the C7 5 D(C7 5 )~ 1 = D* symmetry of the Dirac 
operator is crucial for avoiding the sign problem, but de- 
formations generically break it. There are a variety of 
deformations that prevent baryonic pion condensation, 
and several ways of introducing auxiliary fields. A de- 
formation which turns out to have a simple effect on the 
deformed theory [H] and can be implemented without a 
sign problem in the chiral limit is 



(s^ ab s ab 



pf ab p 
' J n. 



(10) 



Aqcd 

where P a h = qJCqi,. Using Fierz identities, this can be 
written as 



A 



QCD 



(11) 



where there is an implied sum over the color labels 
We then introduce real auxiliary fields that couple to the 
flavor-singlet bilinears qi ■ ■ ■ qj. This allows us to main- 
tain a CDC^ 1 = —D* symmetry for any c, fiB so long 
as m q — 0, avoiding the sign problem. While in practical 
lattice calculations m q > 0, the lack of a sign problem at 
vriq = implies that the phase-quenching approximation 
must become increasingly accurate as m q — > in this 
theory. 

The sign-free deformation may seem peculiar, espe- 
cially since Cd is not positive definite. However, one can 
get a nonperturbative understanding of its effects using 
low-energy effective field theory. The result of this analy- 
sis, which will be presented elsewhere |T3], is very simple: 
the deformation raises the mass of the baryonic pions, 
pushing their condensation point past \ib = sav- 
ing the equivalence. Once [1b ~ ^qcd, baryonic mesons 
with masses ~ A.qcd might condense. Whether this hap- 
pens depends on which states in the deformed theory 
have the smallest mass per U(1)b charge, and will have 
to be resolved by lattice simulations. It would be very 
interesting if the equivalence works through the nuclear 
matter transition. 

More applications. The arguments so far also hold at 
finite temperature. In particular, the details of the chi- 
ral transition can be studied as long as the baryonic pion 
does not condense. So even without the deformation, one 
can gain valuable insights into hot QCD. Our framework 
also gives insights into the behavior of phase-quenched 
simulations. When Nf is even, one can do a projec- 
tion of the SO theory by using a Z4 subgroup of the 
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U(l)i 3 € SU(Nf)v "isospin" flavor symmetry instead of 
U(1)b- The daughter theory is again large N c QCD, but 
now hb in the parent is mapped to an isospin chemical 
potential fii in QCD. The resulting daughter-daughter 
equivalence between QCD with /is 7^ and fii 7^ holds 
for jU < m-n/2. It has been previously noticed that phase- 
quenching QCD with ixb < seems to be a good 
approximation for some observables, and that the phase- 
quenched theory is just QCD with fXj 7^ [20 . The 
daughter-daughter equivalence guarantees that at large 
N c phase quenching is exact for observables with zero 
baryon and isospin charges, giving additional insights 
into the behavior of this approximation [2T] . The equiv- 
alence can be extended to a holographic setup [22], and 
the coincidence of the phase diagram in the barvonic[2"5] 
and isospin[53] theories can be seen analytically. 

Outlook. We have proposed a way to dodge the sign 
problem in the chiral limit of large N c QCD by working 
with a large-iV c equivalent SO(2N c ) theory. There are 
many directions for future work, some of which were men- 
tioned above. Of these, tests of the proposal on the lat- 
tice and EFT analysis of the IR physics of the deformed 
theory are perhaps the most urgent. Finally, one might 
wonder if orbifold equivalence can allow one to dodge 
sign problems in other systems, for instance in SYM the- 
ories [35]- This would make the Monte Carlo approach 
to the gauge-gravity duality (see, e.g. [55]) much more 
tractable. 
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